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ABSTRACT 


The paper considers splitting a plane elastic wedge- 
shaped solid through the application of a rigid punch. It 
is assumed that the coefficient of friction on the contact 
area is constant, the problem has a plane of symmetry with 
respect to loading and geometry, and the crack lies in the 
plane of symmetry. The problem is formulated in terms of 
a system of integral equations with the contact stress and 
the derivative of the crack surface displacement as the un- 
known functions. The solution is obtained for an internal 
crack and for an edge crack. The results include primarily 
the stress intensity factors at the crack tips, and the 
measure of the stress singularity at the wedge apex, and 
at the end points of the contact area. 
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1. INTRODUCTION 


In structural and mechanical design the importance of 
stress concentration around notches has long been well rec- 
ognized and the related theoretical problems have been ex- 
tensively studied (e.g. [1]). Conforming to a sound design 
practice, in most of these studies it is assumed that the 
radius of curvature of the notch is greater than zero. How- 
ever, in failure studies where normally brittle fracture 
would be the expected mode of failure, it may be convenient 
and perhaps even necessary to study the limiting case of 
the problem in which the radius of curvature of the notch 
is zero. This problem too has been studied rather exten- 
sively by assuming that the relevant portion of the solid 
may be approximated by an infinite elastic wedge with a 
given notch angle (e.g. [2,3], for additional references see 
[4]). Most of these studies deal with the solution of the 
related traction and/or displacement boundary value prob- 
lems. However, in real problems the external load is 
applied to the component usually through another solid con- 
tacting the wedge. In these problems if the contact region 
is sufficiently close to the apex of the wedge, it is clear 
that the stress concentration factor or the parameter rep- 
resenting the strength of the stress singularity would 
depend on the details of the distribution of contact stresses 
as well as on the resultant forces. Such contact problems 
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for elastic wedges without or with taking into account the 
effect of friction were considered in [4] and [5]. A 
slightly different problem dealing with the process of 
wedge-splitting of a semi-infinite strip was considered 
in [6]. 

The results given in [4] and [5] may be used in 
brittle wedge-shaped solids to study the initiation of 
fracture at the wedge apex. To deal with the problem 
after fracture initiation one needs the solution of the 
corresponding crack-contact problem. For the symmetric 
geometry shown in Figure 1 and under symmetric loading 
conditions the problem is considered in this paper. In 
formulating the problem it is assumed that the driving 
punch through which the load is applied is rigid and the 
coefficient of friction along the contact area is constant. 
More recently, the traction boundary value problem for an 
infinite wedge containing an edge crack was considered in 
[7], where the problem is solved for uniform pressure on 
the wedge and the crack surfaces. The nature of the ex- 
ternal loads considered in [7] however, appears to limit 
the practical usefulness of its results. 

2. FORMULATION 

Consider the plane elastostatic problem for the in- 
finite wedge-shaped domain of angle 2 0 Q shown in Figure 1. 
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The external load P Q is applied to the medium through a 
rigid punch with known profile. It is assumed that 
(6*0, 8=ir) is a plane of symmetry with respect to geometry 
as well as loading. If the magnitude of the load P Q is 
increased beyond a certain value, then a crack may develop 
along 8=0 plane. This may be the consequence of cyclic 
loading or, in "brittle" materials such as rocks, ceramics, 
and a rather large variety of cast metal alloys, it may 
result from the static application of the load. Generally, 
the crack would start from the apex of the wedge as an 
"edge crack." However, the problem will be formulated for 
the internal crack geometry shown in Figure 1 and the edge 
crack will be treated as the limiting case for c+0. In 
this problem, because of symmetry it is sufficient to con- 
sider one-half of the medium only. First the problem will 
be formulated for the general case in which the normal and 
tangential components of the unknown contact stress vector 
are independent of each other. This formulation is neces- 
sary to solve the problem in which due to a very high co- 
efficient of friction, no sliding may take place between 
the punch and the elastic wedge. In this case the plane 
elasticity problem must be solved under the following 


boundary conditions: 




a 08^ r ’ 0 O^ = 0 ’ 

0<r<a 

, b<r<«> , 

(la) 

3r 4y 

9] (r) 

, a<r<b , 

(lb) 
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a re ( r * 9 o ^ = 0 » 0<r<a 


b<r<°° 


(2a) 


) 


u r (r,e Q ) = ^ g 2 (r) , a<r<b , (2b) 

o re (r,0) * 0 » 0<r<» , (3) 

u 0 (r,O) = 0 , 0<r<c , d<r<°° , (4a) 

o q £ r,0) = g 3 (r) , c<r<d , (4b) 


where y is the shear modulus, <=3-4v for plane strain and 
k=(3-v)/(1+v) for generalized plane stress, and g^ , g 2 , and 
g 3 are known functions. 

The problem is formulated by using the Mellin trans- 
forms. Defining the stress and displacement combinations 

a(r,9) * a r0 + ia QQ , v(r,0) = ~ (u r +iu Q ) , (5a, b) 

for an elastic wedge one may easily obtain (see, for example, 

[ 8 ]) 

M[r z o] * 2i ( s+1 ) [Ase 1 se +■ B(s+1 Je 1 < s+2 > 6 

- §e- 1(s+2)e ] . (6) 

M[r 2 v] = . Ill [Ase 1 ' s9 ♦ B(s*l)e i(s+2)9 

t KBe- i(5+2)9 ] . (7) 

where for a given function f(r) the Mellin transforms are 
defined by 
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f(r) 


1 

2 7T i 


• c 0 + i °° 
C Q ~ 1 00 


M[f]r~ s ds , 

( 8 ) 


provided the strip of regularity containing the constant c Q 
is selected in such a way that, considered together with the 
behavior of f(r) as r-*- 0 and r-*=°, the integrals exist. In 
(6) and (7) the complex functions A ( s ) and B(s) are unknown 
and are determined from the four boundary conditions speci- 
fied along 0=0 and 9 = 6 0 . In the problem under consideration 

there are one homogeneous (i.e., condition (3)) and three 
mixed boundary conditions (i.e., (1), (2) and (4)). Thus, 
substituting from (6) and (7) into (l)-(4), one could elimi- 
nate one of the four real unknown functions and obtain a 
system of three simultaneous dual integral equations for the 
remaining three unknowns. Since the problem involves singu- 
larities with powers other than + 1/2, this would not be the 
proper way to try to solve the problem. On the other hand, 
defining 


f 1 ( r ) = a 69 ( r > 9 0 ) > 0<r<°° , 

(9) 

f 2 ( r ) = a r9 ^ r ’ 9 o ^ » 0<r<°° , 

(10) 

f 3 (r) = TT7 TF u e< r -°> • 0<r< " - 

(ID 


and replacing the conditions (1), (2), and (4), respectively 
by (9), (10), and (11), one may easily determine the complex 
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functions A ( s ) and B(s) and, as a result, all the relevant 
field quantities such as a(r,9) and v(r,9) in terms of f-^ , 
f 2> and f 3> From (la), (2a), and (4a) we observe that f^ , 
f 2 , and f^ are zero on the infinite portion of the interval 
0<r<o°, and are unknown on the remaining part. Using now 
the boundary conditions (lb), (2b), and (4b) which have not 
yet been satisfied, we could obtain a system of integral 
equations having finite supports to determine these func- 
tions. The analysis leading to the integral equations is 
somewhat lengthy and will not be given in this paper. The 
derivations and the asymptotic analysis follow very closely 
and partially repeat the procedure outlined in [4], [5], 
and [6]. Thus, for the problem described in Figure 1 the 
boundary conditions (lb), (2b), and (4b) give the following 
integral equations to determine the unknown functions f -j , 
f 2 , and f^: 

g,(r) = -<£jf) f 2 (r) + 1 ♦ k n (r,t)]f,(t) 

1 1 f d 

+ k 12 (r,t)f 2 (t)j>dt+ - k 13 (r,t)f 3 (t)dt , 

a<r<b , (12) 

g 2^ r ^ = (f+f) M r ) + * J a < [ k 21 Ct) 

+ [~ + k 22 (r,t)]f 2 ( tj^dt 
, f d 

+ - k,-(r , t)f,( t)dt , acrcb , (13) 

* Jc 


- 6 - 



9 3 < r > 


1 


IT 


r b 

C k 31 (r » t) f i ( t) + k 32 (r,t)f 2 (t)]dt 


+ 


1 

TT 


rd 


[-L + 

L rp 


k 33 (r,t)]f 3 ( t)dt 


c<r<d , 


where 


P * log(t/r) 


and the Fredholm kernels k . , ( i , j = 1 , 2 ,3) are given by 

* J 


its i n 2 9. 


k ll (r,t) = " r( 29 +s in20 ) + 
'0 0 ' 


00 cosh29 o y-cos20 o 

„ ( DT71 


si npy 


dy 


k ]2 (r,t) = 


irsi n0 rt cos0„ 
o o 

r(20 o +s1n20 o ) 

« sinh20 n y 


r°° sin20 Q 
rD (y ) 


s i npy dy 


(- 


UT T) 


-1) ^dy , 


k 13 (r,t) = 2 


« cos0 o si nh0 o y+ys i n6 o cosh0 o y 

FDT71 


cospy 


*<21 ( r » t) = ki 9 ( r , t ) • > 


12 ' 


ircos z 9 c w 

k 22 (r,t) = ’r(20 o +sin20 o I + J 0 ( 


cosh29 o y+cos20 o 


TOT 


-l 


X^dy , 


k 23 (r,t) = 2 


» s i n9„ s i nh0 v 
o o 

FoTyJ 


(cospy - ysinpy)dy 


k 31 ( r » ^) ” k n , ( r , t) 5 


13 


(14) 


(15) 


- 1 ) 


dy , 
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k 32 ( r » t ) * 2 


c°° sin0 Q sinh0 o y 
rO(y ) 


(cospy + ysi npy )dy , 


k 33 (r,t) 


f - coshze y-l-y'O-cosZe ) $1 

j ' itttj ’> — r - d * • 


( 1 6a-i ) 


D(y) = sinh20 Q y + ysin20 Q 


(17) 


From the definition of f 3 as given by (11) and the 


boundary condition (4a) it is clear that for c>0 f 3 

must 

satisfy the following si ngl e-val uedness condition: 



■d 

f,(t)dt = 0 
c J 

(18) 

A1 so , 

the force equilibrium along the contact area 

requires 

that 




■b 


« 

f, ( t) d t = -P , 
a 

(19) 


r b 



a f 2 ( t ) d t = -Q , 

(20) 


where P and Q are the normal and the tangential components 
of the resultant contact force. Thus, the integral equa- 
tions (12-14) must be solved under conditions (18-20). 

In the case of constant coefficient friction n 

Q = nP , (21) 

P = P Q /[2si n(7r-0 o ) + 2qcos(7r-e o )] , 


( 22 ) 



where P Q is the driving force shown in Figure 1. In the no- 
slip case (i.e., for perfect adhesion) fj and and conse- 
quently, P and Q are independent and are determined from the 

following equilibrium and kinematic compatibility condi- 
* 

tions : 

Psin(7r-e 0 ) + Q co s ( 7T- 0 Q ) = P Q /2 , (23) 

tan(ir-0 Q ) = u Q (a,e 0 )/u r (a ,e Q ) . (24) 

Examining the behavior of the term 1/rp in the kernels 
i t may be shown that 

rp rl og ( t/r ) r ^t _ 1 j j| n + 1 ^r ^ 

= trp [1 + 0(1 -1)] . (25) 

Thus, the integral equations (12-14) are singular with Cauchy 
type kernels. Also, one may note that the dominant part of 
(12) and (13) is identical to the coupled integral equa tions 
for an elastic half plane loaded by a perfectly adhering rigid 
punch, and that of (14) is the integral equation for the crack 
in an infinite plane. 

If the contact under the punch is frictionless, then 

* 

Note that in this case the contact surface moves by a rigid 
body displacement parallel to the 0=0 plane. Therefore, the 
kinematic relation can be written at any point on the con- 
tact surface. 
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f 2 ( ) = 0 and (12) and (14) give the system of singular in- 
tegral equations to determine the unknown functions f^ and 
f^. In this case the additional conditions are (18) and 
the following eqiulibrium relation: 

r b 

f, ( t)dt = -P /2sin( 1T-0 ) . (26) 

J a 

On the other hand, if the punch is driven with constant co- 
efficient of friction n» we have 


f 2 ( t) = ( t) 


g 1 (r) = — y f -j (r) + - 


rb 


^rp + k ll ( r * t ) + nk 


+ 1 

TT 


rd 


3 ( r,t)f^( t)dt , a<r<b 


9 3 < r > * ? 


[ k 3 1 ( r , t ) + n k 32 (' r » t ) 1 f \ ( t)dt 


+ 1 

IT 


, rd 


T 33 

Y = n ( k- 1 )/(<+!) 


[rrr + k,,(r,t)]f,(t)dt , 


(27) 

1 2 ( r , t ) ]f 1 ( t)dt 

(28) 

c<r<d , (29) 
(30) 


where the kernels k. • are given by (16). In this problem 

* 0 

(28) and (29) must be solved under conditions (18) and (22). 
3. STRESS SINGULARITIES 


Referring to Figure 1 the stress state is expected to 
be singular at the crack tips c and d and at the wedge apex 0. 
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Defining the stress intensity factors at the crack tips by 


k(c) = 


1 im 
r-*-c 


/2(c-r) o 00 (r,O) 


k ( d ) = 1 im /2( r-dj a (r, 0 ) 
r-*d 00 


(31a, b) 


and noting that in (14) g 3 (r)=a 00 (r ,0) outside as well as 
inside the interval c<r<d, one may easily show that only 
the dominant term in the integral equation would contribute 
to the singularity giving [9] 


k ( c ) = 1 i m / 2 (r-c) f ( r ) , (32) 

r+c 

k(d) = -lim /2 ( d -r ) f^(r) . (33) 

r-*-d 

If c>0, the wedge apex is also a point of stress singularity. 
Perhaps the easiest way to extract this singularity would be 
going back to ( 6 ) and expressing a 0Q in terms of f ^ , f 2 , and 
f 3< Thus, after some simple manipulations we find 

/■b 

ro 0 Q (r, 0 ) = I [K-j ( r , t) f -J ( t) + K 2 ( r , t) f 2 ( t) ]dt 



+ 

rd 

K 3 ( r , t) f 3 ( t) dt , 

(34) 

Kj ( r , t) 

1 

2 7T i 

rc + i°° t s+1 F . ( s ) 

0 . (?) -Arrr ds > 

C Q -1°o 

(35) 

A(s) = 

(s+l)sin29 0 + sin2(s+l)9 0 , 


F^s) = 

( s +2 ) s i n ( s +2 ) 0 Q - ssins9 0 , 
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F 2 (s) = scos(s+2)9 0 - scoss6 0 , 

F 3 (s) = s(s+2) - (s+1 ) 2 cos26 0 + cos2(s + l )0 , (36a-d) 

where the strip of regularity for Mel 1 i n inversion integrals 
is given by [4,5] 

Re(s_ 1 )<Re{s) = c Q <-l (37) 

s_.j being the first root of A(s)=Q to the left of the line 
Re(s)=-1. For small values of r, t>r and consequently to 
evaluate the integrals in (35) the contour must be closed to 
the left. Therefore, from (34) and (35) it is seen that the 
leading term in the infinite series giving the stress will 
be 

a Q0 (r,O) - r _a) , u) = 2 + s_i , -2<s_ 1 <-l . (38) 


Thus, defining the strength of the stress singularity by 


k ( 0) = lim r w 
r+0 


a ee ( r» o) 


(39) 


we obtain 


MO) • 


ATi-f] 


1+s.-] Fp ( s , ) b l+s_i 

f-,(t)t dt + j a f 2 (t)t dt 


. F 3< s -1> f d 

* aTTTT L 


1 +s _] 

f 3 ( t)t dt 


(40) 


where 



A 1 ( s ) = sin20 o + 28 q cos2(s + 1 )8 


(41 ) 


4. SPECIAL CASES: 0 Q = it , 8 Q = tt/2 


In the special cases of the plane with a semi -infi ni te 
crack (0 Q =Tt ) and the half plane ( 8 Q = -rr/ 2 ) , the kernels can 
be evaluated in closed form, simplifying the integral equa- 
tions quite considerably. This requires the evaluation of 
a number of relatively simple Fourier integrals (see for 
example, [5]). Thus, for example, for 0 o = 7r the integral 
equations (28) and (29) become 


g^r) = -yf-j (r) + ( t/r) 1 /2 f 1 ( t)dt 


( t/r) 1 /2 f 3 (t)dt , a<r<b 


2 

TT 


g 3 (r) . - I 


+ 1 

TT 


rd 


a 

rd 


(t/r) 1/2 f,(t)dt 
jTp ( t/r)' ^ 2 f j( t)dt , c<r<d 


Similarly for 0 Q = tt/2 we find 


r b 


g^r) = -yf-j ( r) + \ t*- 7 2 f ]( t ) dt 


+ 1 

tt 


rd 4rt 2 


( t 'i + r 2 ) T f 3^ t ) dt » a <r < b 




♦ i 

TT 


(42) 


(43) 


(44) 


fe + ii fwr i 3f3(t) d t , c<r<d 


(45) 
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Equations (42-45) reduce to the integral equationsfor 
the well-known crack and contact problems if the crack and 
contact regions are taken sufficiently far apart. For ex- 
ample, in (42) and (43) if c -*», d-*» with d-c = finite, the 
coupling terms vanish and the uncoupled equations become 
that of a punch in a semi -i nf i ni te crack given in [5] and 
that of a crack in an infinite solid (which has only a 
simple Cauchy kernel). Similarly, in (44) and (45) for 
c-*», d-*», d-c = finite the coupling terms and the second 
term in (45) vanish giving again the known equations of 
two symmetric punches on a half plane and a crack in an 
infinite plane. Also, for a-*», b-*», b-a = finite, the 
coupling terms again vanish, (44) reduces to the integral 
equation for a single punch on a half-plane (which would 
now have in addition to the term yf^ only a Cauchy kernel), 
and (45) becomes the integral equation of a semi -i nf i ni te 
plane having a crack perpendicular to the boundary (see, 
for exampl e [10]). 

In a half plane subjected to compression through rigid 
punches, since normally 9=0 plane would be under compression, 
the case of 6 Q = tt/ 2 may not be very practical. 
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5. THE NUMERICAL SOLUTION 

The dominant part of the singular integral equations 
(14) and (29) is of the first kind, is not coupled with that 
of the remaining equations in the system, and has only a 
simple Cauchy-type kernel. Since the singular behavior of 
the solution is dependent on the dominant part of the inte- 
gral equations only, the unknown function f.j(r) would have 
as expected, standard square root singularity. The dominant 
parts of (12) and (13) are coupled and the integral equations 
are of the second kind. By defining 

f (r) = f-|(r) + 1f 2 (r) , a<r<b (46) 

the dominant parts can be combined and the solution may be 
expressed as 

f(r) = F(r) (b-r) a (r-a) 0 , -l<Re(ct,S)<l . (47) 

where F(r) is bounded in a^r^b. 

Similarly, from (29) it follows that [9] 

f«l (r) = F ] (r)(b-r) a (r-a) 0 , -l<Re(a,8)<l . (48) 

where, again F-|(r) is bounded in a^r^b . The constants a 
and 6 may be determined following the standard function- 
theoretic method (see, for example, [9]). The numerical 
values of a and B and the index of the integral equations 
(12), (13), and (28) would depend on the conditions of con- 
tact and the profile of the punch. These values are tabu- 
lated in Table 1 where n = » corresponds to the perfect 
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adhesion case (in which f-j and are independent) and the 
constants Q and (jj q are defined by 


tana = v 

n(<- i ) 


O). 


_ log < 


“o 2-rr 

One may note that in the case of constant coefficient of 
friction a and 8 are always real, whereas in the perfect 
adhesion case they are complex and hence the singularity is 
of oscillating type. 


The numerical results given in this paper are for the 
case of constant coefficient of friction. The corresponding 


Table 1. Values of a and 8 




■1 

= 0 

n = constant 

n = 

OO 


Index 

a 

. 6. 

a 

8 

a 

B 



1 

1 

"2 

1 

"2 

B 

_a 

~ IT 

RBB 

"F +iu o 

J 


-1 

1 

7 

1 

2 

a 

IT 

1 Q 
t r 

2~ la) o 

EH 

J 


0 

1 

"I 

1 

2 

B 

i _a_ 
TT 

m 

b u o 

j 


0 

1 

2 

1 

‘2 

n 

IT 

_a 

r 1u o 



integral equations are (28) and (29) which must be solved 
under conditions (18) and (19). The solution is obtained 
by using the Gauss -Chebyshev and Gauss-Jacobi integration 
formulas described, for. example, in [9]. Before using 
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these formulas the intervals in (28) and (29) are normalized 
by defining 


Thus 

by 


_ 2r 


b+a 


2t 


b+a 


A 1 b-a ' b-a ’ *1 b-a " b-a ’ a< ( r » t ) <b 

Y _ 2r d+c _ _ 2t d+c , 

x 2 - » s 2 ' ‘ ’ c<(r, t)<d 

♦l ( s 1 ) = f 1 (t)/p Q , <j> 2 (s 2 ) = f 3 ( t ) /P Q , 

( x 1 ) = g 1 (r)/p Q , ^ 2 (x 2 ) = 9 3 (r)/P 0 , 

p 0 « 2P/ (b-a) , -l<(x ] ,x 2 ,s 1 ,s 2 )<l . 


(50) 


, equations (28), (29), (19) and (18) may be replaced 


1 

^l( x l) = - Y<f>i ( *i ) + ^ ^1 1 ^ x i »s-| ) 4 >i ( s-j ) d s i 


-1 


+1 

TT 




^ 1 2 ^ x 1 ,s 2^2^ s 2^ s 2 * ~ ^ < x -| < 1 


> 


^ 2 (x 2 ) 



rl 




(x 2 , 

( x 2 » 


S 1 

s 2 


)<j» 1 (s 1 )ds 1 
)4> 2 (s 2 )ds 2 


1 

<^1 ( s l ) d s -| = -1 » 


rl 


- 1 < x 2 < 1 


J ^2 ^ s 2 ^ S 2 = ^ • 

Where the kernels h... (i,j=l,2) are obtained by 

* v 

the respective terms in (28) and (51a), and (29) 
Defining now 


, (51 a ,b) 


(52 a ,b) 

compari ng 
and (51b). 


«t> 1 (s 1 ) = G-j (s ] )(l-s 1 ) a (l+s 1 ) 6 , 

I 

<f> 2 ( s 2 ) = G 2 ( s 2 )(l-s|)“2 , (53 a,b) 
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and using the quadrature formulas given in [9], (51) and 
(52) become 


k^ H k h n^11 » s lk^ G l t s lk J + W k h 12^1i > s 2k> G 2 (s 2k^ 

= ^(Xj.) , 

k S 1 C H k h 2 1 ^ x 2i ,s lk^ G l ^ s lk^ + W k h 22^ x 2i * s 2k^ G 2^ s 2k^ 

= 7r '^2 ^ x 2i ^ » * = 1 n - 1 , (54 a,b) 

k-1 H k G l ^ s 1 k ^ = _1 ’k = l W k G 2^ s 2k^ 0 » ^ 55 a,b ) 

where it is assumed that the punch has "sharp" corners. 

Hence, a = (£2/7 r) - 1 and 3 = -£2/7r(see Table 1). The weights, 
integration points, and collocation points which appear in 
(54) and (55) are given by [9] 



X 


2n -2 + ct+B 

p ( l+o, 1 + 3) / . \ d p (l+-ct,l+3) / _ \ 

K n-3 ' S 1 k 'ds^n-2 is lk ; 

i 

k = 2 , . . , n - 1 ; 


= 1 s = -1 p ( 1 +a , 1 +3 ) / _ \ _ g 

11 1 ’ s ln 1 * P n-2 ls lk ; U 

k = 2 , . . , n - 1 ; 
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pM-a.-’-e)^) . p ( j = ] „ - 1 ; 

W 1 = M n = 2(n-1 ) ■’ W k = fTT ■ k “ 2 1 ; 

k-1 

s 2 — cos (ff~ ) » k - 1 , . . , n > 

*2 ^ “ cos^tt 2 p “ 2 ) » i - 1 . • • » n - 1 . (56) 

Equations (54) and (55) give 2n linear algebraic equations 
to determine the unknowns S^(s^ k ) and G 2 ^ s 2k^ > k = 1 

> • • ) n • 


In all the examples given in this paper it is assumed 
that the crack surfaces are traction-free. Therefore, 
g 3 = 0 or ~ °* Also, in a flat-ended punch with 

sharp corners ^(x-j) too is zero. 


For a blunted punch which consists of a flat face and 
a rounded nose with radius R in the contact region the in- 
put function g-j(r) is given by (see the insert in Figure 7) 


g(r) 


4 yi 

1+K 


r - a o 

R 


a< r <a o 


(57) 


10 , a Q <r<b 

where a Q , b, and R are known constants and a is unknown. 

In this case the index of the integral equation (28) is 
zero (see Table 1), meaning that no additional condition is 
required for a unique solution. The unknown constant a is 
then determined from the equilibrium condition (19). To 
solve the problem, first the equations are again normalized 
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by defining the new variables x-| , s^, x 2 , s 2 as given by 
(50) and letting 

4>-, Cs 1 ) = f 1 (t ) , d > 2 Cs 2 ) = f 3 Ct ) , 'I'-j (x ] ) = g-| (r) , 

i>2 ( x 2 ^ = 9 3 C r ) . (58) 

Then, the system of 2n algebraic equations consisting of 
(54) and (55b) and 

G-, Cs ln ) = G-j (-1 ) = 0 (59) 

is solved for G-j(s lk ) and G 2 (s 2k ) for a specified value of 
the constant a. The corresponding value of the resultant 
force P is then obtained (in an inverse manner) by using 
the equilibrium condition (19) or 



In the case of the edge crack, i.e., for c=0, basically 
the integral equations (12-14) or (28,29) and the kernels 
(16) remain the same. However, in this problem even though 
the equilibrium conditions such as (19,20) and (26) must still 
be satisfied, the single-valuedness condition (18) is no 
longer valid. The reason for this, of course, is that at 
r =0 u Q (r, 0 ) is no longer zero and f 3 (r)= 3 u Q / 3 r is nonzero, 
finite, and unknown. Also note that at r=0 the asymptotic 
problem is that of an elastic wedge of angle 0 0 » ( 0 < 6 0 <- ^ ) for 
which j ( r , 8 ) ~p i j+ 0 ( r a ) as r-*-+ 0 , (i,j=r, 8 ; 0 < 8 <e o ) where p^ 
is a finite constant which is zero for wedges with stress- 
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free boundaries and the real constant a is positive.* In this 
case the numerical solution of the problem can be obtained by 
following the procedure outlined in this section with a minor 
change, namely by excluding the s i ngl e-val uedness condition 
(55b) and by replacing it by the condition S 2^ s 2n^ =0, Thus » 
considering (53b), at s 2 = -l (or at r = 0) <j> 2 (s 2 ) becomes inde- 
terminate and may be evaluated by extrapolation (if needed). 


6. RESULTS AND DISCUSSION 

The results for the wedge problem without the crack 
were given in [4^ and [5]. However, in order to give some 
Idea about the effect of the wedge angle and of the coef- 
ficient of friction on the stresses which may cause the 
formation of a crack, the stress distribution along 9=0 
plane is calculated. The results are shown in Figure 2. 
Based on these results, one could make two important obser 
vations, namely the existence of friction would inhibit 
crack initiation and, for relatively small wedge angles 
29 . as r increases the cleavage stress c QO (r,0) would 

0 t7 V 


The s tatement in Ref. [7] (p. 619) regarding this point, 
namely "A simple residue calculation shows that the stresses 
are 0(r~ s *'') as r-n-0, where s is the first zero of a(s, 0 q ) 

= ssin29 +sin2s0 to the left of |ero. It is easily shown 
that as 0° increases from 0 to it, s increases from -3/2 to 
-1/2," is°clearly incorrect. It is valid only for a wedge 
of angle 20 o under symmetric loading, that is, in the present 
case, only for the uncracked wedge (see, equations 36a, 37, 
and 38, where s+1 replaces s in [7]). Upon introducing the 
crack, the total wedge angle becomes 0 O , ( 0< 0 O <tt ) and the 
stress singularity disappears. 
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change sign and would become negative. This means that 
for such wedge angles, even tf the crack is formed, it 
would not propagate very far. The calculated stress inten- 
sity factors given in this paper also support this conclu- 
sion. 

In the elastic wedge the practical problem is, of 
course, the edge crack problem. This is strongly demon- 
strated by the results for an internal crack given in 
Tables 2 and 3, and Figure 3. The stress intensity factors 
k(o), k(c), and k(d) given in the tables and in the figure 
are defined by equations (39), (31a) and (31b), respectively. 
As the crack distance c/£ increases, as indicated by Figure 
3, k(o) quickly converges to the values given in [5] for 
the uncracked wedge. In the figure and in the tables k(o) 

Is normalized with respect to P / ( Tre ^“ w ), e = (b+a)/2 and 
w = 2 + s_-j (see equations 38 and 39). The results given 
in tables 2 and 3 and Figure 3 show that as c+o, k(o)-»-«> 
and k(c)-«“. Hence in brittle solids the most likely loca- 
tion of crack initiation would be the wedge apex. The tables 
also show the strength of singularity at r = a and r = b 

for the contact stress a QQ (r,0 ). Recalling the form of the 

0 0 0 

contact stress given by (48), these stress concentration (or 
Intensity) factors are defined by 

k (a ) = - 1 i m /? ( r-a) ' e a 00 ( r , 0 Q ) , 

k(b) = -lim (b-r)' a a 00 (r,0 o ) . (61 a,b) 
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Table 2. The results for the flat punch with sharp corners 

and the internal crack (see the insert in Figure 3), 
* = ( b - a ) / 2 , e = (b+a)/2, m = (d-c)/2, v = 0.3, n = 




CO 

in 


VO 

in 

CO 

in 

r— 

I s ** 

r** 

CO 

in 

ID 

CO 


■ — tc«J 

1— 

cn 

CM 


00 

<n 

co 

in 

CO 

CO 

CM 

CO 

p— 

CM 

T5 

E 

Cn 

in 

cn 

CO 

CM 

in 

o 

CM 

CO 

o 

CO 

r— . 

CO 


" 

t= 

O 

CM 

CM 

co 

r— 

CO 


in 

o 

CM 

CM 

r— 

CM 

CO 


' ** 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 



o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 


CL 

















a> 

cn 

vo 



cn 

vo 


vo 

o 

m 

00 

VO 

CM 


■— |CM 


co 

cn 


r-*. 

r— 

00 


c 

in 

co 


OO 

00 

U 

E 

CM 

CM 

CM 



oo 

co 


r— 

r** 

00 

in 


CM 

'W 

t= 

r— 

CM 

cn 


i— 

VO 

co 


r— 

CO 

o 

r— ■ 

r— 



« • 

• 

• 

« 


• 

• 

• 


• 

• 

• 

• 

| • 

• 



o 

r— 

CM 


o 

r— 

CO 


o 

o 

CM 

o 

i 

CM 


a. 
































3 















s 

+ 

CT> 

r**. 


CM 

cn 


o 

in 

r— 

00 

VO 

00 

f— 

00 

JO 

f— 

VO 

o 

vo 


in 

oo 

00 

in 

I s *. 

1 — 

CO 

in 

CO 

in 



O 

CM 

CM 


CM 


in 

CO 

in 

in 

in 

i — 

CM 

CM 



cn 

cn 

cn 

cn 

00 

CO 

CO 

CO 

cn 

cn 

cn 

cn 

cn 

cn 


*w* 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


**■**„ 

© 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


CL 
































<n 















' 

+ 

CM 

o 

o 

VO 

in 

in 

^1“ 

co 

CO 


r^* 

cn 

CM 

CM 

<T3 

r— 

in 

o 


cn 

o 

vo 

VO 


in 

00 

in 

vo 

VO 

CM 

w - 

oi 

in 

CO 


I s *. 

00 

CO 

r— 

in 

in 



o 

cn 

cn 


t= 


*— 

r— 

o 

CM 

CM 

CM 

r— 

o 

o 

o 

f— 

o 

o 



• 

• 

# 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 


N. 






^* 


r— 

I— 


r— 

r— 

r— 

r— 


CL 
































3 
















1 

r— 




in 

cn 

cn 


r** 

VO 

r— 

cn 


r— 

o 

U— 

cn 

CM 

VO 


CT> 

cn 

LO 


co 

VO 

in 

VO 

cn 

CM 

>— * 

0) 


CO 

CO 


r**. 

CO 

o 




o 

CM 

o 

in 


!= 


cn 



cn 




r*. 

cn 

in 

cn 


vo 



• 

• 

• 


• 

• 



• 

• 

• 

• 

• 

• 



O 

r— 

in 


o 

CM 

VO 


o 

r— 

in 

o 

CM 

VO 


CL 



























r* 







r— 

o 



r— 

o 



t — 

o 


r— 

o 




o 

o 



o 

o 



o 

o 


o 

o 

T3N 

• 

• 

• 


• 

• 

• 


• 

• 

• 

• 

• 

• 



CM 

r— 



CM 




CM 

r_ 

r_ 

CM 






















r— 

o 



l— 

o 




o 


r— 

o 




o 

o 



o 

o 



o 

o 


o 

o 

<j|<* 

• 

• 



• 

• 

• 


• 

• 

• 

• 

• 

• 



r* 

o 

o 

o 

r— 

o 

o 

o 

' 

o 

o 

r— 

o 

o 

X3M 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

in 

in 

in 

in 

in 

in 



— * 


r- 

*— ■ 



r- 

co 

CO 

oo 

co 

co 

co 



0 

0 

0 

o 

o 

0 

0 

o 

o 

o 

o 

o 

O 

o 


o 

o 

o 

o 

o 

in 

in 

LO 

in 

o 

o 

o 

in 

in 

in 

® ! 

in 

in 

in 

in 

vo 

VO 

lO 

vo 

in 

in 

in 

vo 

vo 

in 

_____ 

1 

p— 


r- 

• — 

» 


* — 


• — 

r— 


1 

' 

»— 


- 23 - 




Table 3. The results for the flat punch with sharp corners 

and the internal crack (see the insert in Figure 3), 
l = ( b - a ) / 2 , e = (b+a)/2, m = (d-c)/2, v = 0.3, n = 
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Without covering a wide range, the tables show the trends of 
various stress concentrations for varying crack distance 
(d+c)/2, punch distance (b+a)/2, wedge angle 20 o> and coef- 
ficient of friction n« 

The results obtained for the edge crack (c=o) are 
given by Figures 4-7 and Tables 4-9. For fixed relative 
dimensions d = a = i , (b-a=2M and for various values of 
the coefficient of friction n» Figure 4 shows the effect of 
the wedge angle 20 Q on the stress Intensity factor k(d). 

It is seen that k(d) decreases as n increases or as 0 Q 
decreases, becoming eventual ly zero on a critical line in 
n vs. 0 Q plane. For fixed values of wedge angle (0 o = 15O°), 
crack length (d/a=l ) , and coefficient of friction (n=0, 
n=0.5) , the effect of the contact area as measured by b/a 
is shown in Figure 5. As b increases there is a slight 
reduction in the stress intensity factor k(d). Figure 6 
shows the effect of the coefficient of friction n on the 
stress intensity factor k ( d ) for fixed relative dimensions 
and for 0 Q = 150°, 165°, 180°. For 0 Q = 180° the change in 
k(d) (which is quite insignificant) comes primarily from 
the effect of n on the contact stress cr QQ (r,0„). 

For the flat punch with sharp corners more detailed 
results are given in Tables 4-7. The negative k(d) values 
found for 0 Q = 120° in Table 4 indicate (from the facture 
viewpoint) the beneficial effect of friction and should be 
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Table 4. The results for the wedge containing 
an edge crack loaded by a flat punch 
with sharp corners, v = 0.3, a/£ = 1, 
b/A = 3, A = ( b - a ) / 2 . 


n 

e o 


■jiBSI 



EZQii 

■Hal 

0 

120° 

0.8832 

1 .0999 

0.1133 

0 

150° 

1.1267 

0.9058 

0.6876 

0 

165° 

1 .1826 

0.8600 

0.8123 

0 

172.5° 

1.1946 

0.8501 

0.8391 

0 

180° 

1.1981 

0.8472 

0.8472 

0.2 

120° 

WKSSM 

1.1526 

-0.1703 

0.2 

150° 

1.0796 

0.9442 

0.5439 

CM 

• 

O 

165° 

1 .1573 

0.8855 

0.7426 

0.2 

172.5° 

1.1780 

0.8699 

0.8030 

0.2 

180° 

1 .1887 

0.8620 

0.8419 

in 

o 

120° 

0.6510 

1.2197 

-0.5834 

o 

• 

cn 

150° 

1 .0003 

0.9935 

0.3323 

0.5 

165° 

1.1115 

0.91 74 

0.6396 

0.5 

172.5° 

1 .1459 

0.8938 

0.7495 

0.5 

180° 

1.1676 

0.8790 

0.8340 
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Table 5. The results for the wedge containing 
an edge crack loaded by a flat punch 
with sharp corners, v = 0.3, a/£ = 1, 
b/£ = 3, 6 = 150°, and n = 0. 


d 

k(d) 

k ( a ) 

k ( b ) 

I 

P/l T/t 

P/(tt£ 1+B ) 

p/(TT-e Uot ) 

0.001 

0.8698 

1 .2192 

0.8575 

0.1 

0.8803 

1 .2030 

0.8648 

0.2 

0.8522 

1 .1892 

0.8714 

0.5 

0.7767 

1 .1585 

0.8874 

1 . 

0.6876 

1 .1267 

0.9058 

2. 

0.5803 

1 .0912 

0.9279 

3. 

0.5145 

1 .0708 

0.9416 

4. 

0.4683 

1 .0572 

0.9511 

5. 

0.4332 

1 .0475 

0.9583 

6. 

0.4053 

1 .0403 

0.9638 

7. 

0.3824 

1 .0346 

0.9683 

8. 

0.3630 

1 .0300 

0.9719 

9. 

0.3464 

1 .0263 

0.9750 

10. 

0.3319 

1 .0232 

0.9776 

























































































d 

l 



0.3323 

0.2930 

0.2697 

0.2525 

0.2387 

0.2270 

0.2170 

0.2082 

0.2004 

0.1934 


0.4430 


0.4336 


0.4134 


0.3712 


0.001 

0.1 

0.2 

0.5 


k 


k 


lb)_ 

o 1 +CS 



■msm 

1 .0484 

0.9713 

1 .0397 

0.9748 

1 .0324 

0.9779 

1 .0166 

0.9852 

1 .0003 

0.9936 

0.9817 

1 .0039 

0.9705 

1 .0105 

0.9628 

1.0153 

0.9572 

1 .0189 

0.9530 

1.0218 

0.9496 

1 .0242 

0.9469 

1 .0261 

0.9446 

1 .0278 

0.9427 

1 .0292 




























































disregarded. Tables 5-7 show the stress intensity factors 
as a function of the crack length d. The important conclu- 
sion which may be drawn from these results is that for a 
fixed driving force P Q , k(d) is a decreasing function of d. 
Therefore, the related brittle fracture process is expected 
to be stable. 

Finally, Table 8 and 9 and Figure 7 show some results 
for an elastic wedge containing an edge crack and loaded 
by a blunted punch. The punch profile is rounded in the 
neighborhood of r = a and has a radius of curvature R 
(see equation (57)). The trailing edge of the punch, r = b 
is sharp. Figure 7 shows a sample result giving the resul- 
tant contact force P and the stress intensity factor k(d) 
as a function of the relative change in the contact area 
(a Q -a)/R. The numerical results for two punch sizes are 
given in Tables 8 and 9. 

In calculating the results the basic computer program 
was run by using the input parameters for known special 
cases for the purpose of verification. The comparison, 
whenever possible, was quite satisfactory. For example, 
in the case of an edge crack with 0 Q = 180° the problem 
is that of a wedge-loaded semi -i nfi ni te crack. In the 
punch-loaded wedge problem letting a /£ = 99, b/Z = 101, 
d fZ - 100, Z = (b-a)/2, it was found that k(d)/[P/(ir/cT)] 

= 0.999995 whereas the exact value is 1 (more specifically. 
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Table 8. The results for an elastic wedge 

containing an edge crack and loaded 
by a blunted punch, v = 0.3, n = 0.5 
8 = 165°, a = R, d = R, and b = 2R 


v a 

1+K P 

_ k(b) 

k(d) 

R 

4y R 

P/R a+1 

1 

P / d 2 

0.5 

0.5990 

0.3742 

0.1891 

0.4 

0.4225 

0.3833 

0.1881 

0.3 

0.2710 

0.3932 

0.1871 

0.2 

0.1460 

0.4037 

0.1861 

0.1 

0.0515 

0.4149 

0.1852 

0.05 

0.0186 

0.4210 

0.1847 

0.02 

0.0048 

0.4246 

0.1845 
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Table 9. The results for an elastic wedge 

containing an edge crack and loaded 
by a blunted punch, v = 0.3, n = 0.5 
0 = 165°, a = R, d = R, and b = 3R 


v a 

l+< P 

1 

k(d) 

R 

4p R 

P/R a+ ' 

P/d2 

0.5 

0.7108 

0.271 1 

0.2123 

0.4 

0.5134 

0.2757 

0.2102 

0.3 

0.3366 

0.2802 

0.2084 

0.2 

0.1851 

0.2848 

0.2067 

0.1 

0.0764 

0.2897 

0.2050 

0.05 

0.0246 

0.2922 

0.2042 

0.02 

0.0066 

0.2938 

0.2037 
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k = /2P/(tt/F), c being the distance of the 
wedge force P from the crack tip). 


concentrated 
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Figure 1. Basic geometry of the crack-contact problem. 
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2.0 



Figure 2. Distribution of the cleavage stress in a punch- 

loaded elastic wedge, a /£ = 1, b/l =3,1= (b-a)/2, 
n: the coefficient of friction. 




Figure 3. Variation of the strength of stress singularity 
at the wedge apex with the crack distance, a /•£ 
b /l = 3, (d-c )/l = 1, the normalization factor: 

k = P/ (ire 1 " u ) , l = (b-a )/2 , e = (b+a)/2. 


c\j 





Figure 5. Variation of the stress intensity factor k(d) with 
the punch size, 0 = 150°, d/a = 1, n: the coef- 

ficient of friction. 







k(d ) 


Stress intensity factor k(d) and the resultant 
contact force P as functions of the relative 
change in the contact area. P = 0 for a = a Q , 

is the point of tangency on the punch profile, 
v = 0.3, n = 0.5, 0 Q = 165°, a_ = R, d = R, b 




